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Abstract. SgrA* is the supermassive black hole candidate at the center of the Galaxy and 
an ideal laboratory to test general relativity. Following previous work by other authors, we 
use the Polish doughnut model to describe an optically thin and constant angular momentum 
ion torus in hydrodynamical equilibrium and model the accretion structure around SgrA*. 
The radiation mechanisms are bremsstrahlung, synchrotron emission, and inverse Compton 
scattering. We compute the spectrum as seen by a distant observer in Kerr and non-Kerr 
spacetimes and we study how an accurate measurement can constrain possible deviations 
form the Kerr solution. As in the case of emission from a thin accretion disk, we find a 
substantial degeneracy between the determination of the spin and of possible deviations 
from the Kerr geometry, even when the parameters of the ion torus are fixed. This means 
that this technique cannot independently test the nature of SgrA* even in the presence of 
good data and with the systematics under control. However, it might do it in combination 
with other measurements (black hole shadow, radio pulsar, etc.). 
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1 Introduction 

SgrA* is the radio source associated to the supermassive black hole (BH) candidate at the 
center of our Galaxy [1]. Its mass is around 4 million Solar masses and its distance from us 
is about 8 kpc. This makes SgrA* a peculiar source and an ideal laboratory to test general 
relativity. It is much closer than any other supermassive BH candidate and much heavier 
than Galactic stellar-mass BH candidates in X-ray binaries. It can thus be studied with 
observations not feasible in other cases. For instance, its apparent angular size on the sky is 
expected to be about 50 ^as and we may be able to image its shadow in the near future [2, 3]. 

Astrophysical BH candidates are so called because they can be naturally interpreted 
as the Kerr BHs predicted by general relativity, but direct observational evidence is still 
lacking [4-6]. Today we have robust measurements of the masses of these objects. This 
is enough to conclude that stellar-mass BH candidates in X-ray binaries are too heavy to 
be neutron stars [7] and that supermassive BH candidates in galactic nuclei are too heavy, 
compact, and old to be clusters of neutron stars [8]. After that, we have to rely on the 
validity of general relativity, even if the theory has only been tested in weak gravitational 
fields, mainly with Solar System experiments and radio pulsar data. In general relativity and 
under a set of reasonable assumptions about the matter content and the spacetime structure, 
the final product of a gravitational collapse is a BH [9, 10]. The only uncharged BH solution 
in 4D is the Kerr metric and therefore BH candidates must be Kerr BHs if general relativity 
is correct and they may not be Kerr BHs only in presence of new physics. 

In the past few years, there have been significant efforts to figure out how to test the 
nature of BH candidates in order to confirm the Kerr paradigm. The study of the thermal 
spectrum of geometrically thin disks and the analysis of the iron Ka line are the two main 
techniques capable of probing the spacetime geometry around these objects and test the Kerr 
metric [11-14]. However, the thermal spectrum of a thin disk has a very simple shape and 
there is a degeneracy between the estimate of the spin and of possible deviations from the Kerr 
background. Eventually, even large deviations from the Kerr solutions are allowed [15, 16]. 
The profile of the iron line has a more complicated structure and it is potentially a more 
powerful tool to test the Kerr metric [17, 18]. However, high quality data would be necessary. 
With the available observations, we can conclude that BH candidates cannot be some types 
of compact dark stars or wormholes [19, 20], because the iron lines of these objects would be 
qualitatively different from those observed, but we still cannot exclude large deviations from 
the Kerr geometry [18]. 
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Since a number of studies has shown that the main obstacle to test the Kerr paradigm 
is the strong correlation between the spin and possible non-Kerr features, it sounds natural 
to try to break the degeneracy with a combination of different measurements of the same 
source [21-25]. With this spirit, SgrA* may soon become one of the best objects to test the 
Kerr metric. While there are currently no measurements to probe the spacetime close to 
this object, there are strong motivations to expect that near future facilities will be able to 
provide a number of different data. The detection of the boundary of the BH shadow can 
provide a measurement of the apparent photon capture radius [26, 27], which only depends 
on the metric of the spacetime [28-38] . The observation of blobs of plasma orbiting near the 
innermost stable circular orbit can measure the orbital frequency and effects related to strong 
gravitational lensing [39-42]. The observation of a radio pulsar in a close orbit, but also of 
ordinary stars, can provide a measurement of the spin of SgrA* [43-45]. The combination 
of different measurements seem to be the right strategy to break the degeneracy between an 
estimate of the spin and possible non-Kerr features, and SgrA* may be the best source to 
get several independent observations to combine together [46, 47]. 

In this paper, we want to explore one more approach to test the Kerr metric with 
SgrA*. The accretion structure around SgrA* seems to be a radiatively inefficient advection 
dominated accretion flow [48]. As proposed in Refs. [49, 50], such an accretion structure may 
be described by an optically thin and constant angular momentum ion torus within the Polish 
doughnut model [51, 52] and a 2-temperature plasma in which the emission mechanisms are 
bremsstrahlung, synchrotron radiation, and inverse Compton scattering. Since the structure 
of the torus is determined by the spacetime geometry, the comparison between theoretical 
predictions with observations can potentially constrain the metric around SgrA*. In this 
explorative work, we compute the spectrum of ion tori in Kerr and non-Kerr backgrounds 
and pay some attention to the correlation among the spin parameter of the central object, 
possible deviations from the Kerr solution, and the fluid angular momentum. These three 
parameters are related to the spacetime metric. However, the spectrum also depends on 5 
other parameters. Here we do not study the correlation among all the parameters in the 
model. 

Even when the parameters of the ion torus are fixed, the spectrum is degenerate with 
respect to the spin and the deformation parameters and it is impossible to test the Kerr metric 
without other data. The situation is quite similar to that of the measurement of the thermal 
spectrum of a thin disk of a stellar-mass BH candidate in a binary. When the fluid angular 
momentum is also a free parameter to be determined by the fit, the constraining power of 
the spectrum is substantially the same. In order to test SgrA*, it is necessary to combine 
different observations. The constraints from the spectrum of the accretion structure might 
be used, for instance, in combination with those from the observation of the BH shadow, 
from the detection of a pulsar/stars in close orbits, and from the observations of blobs of 
plasma orbiting near the innermost stable circular orbit. 

The content of the paper is as follows. In Section 2, we review the Polish doughnut model 
in non-Kerr backgrounds. In Section 3, we summarize the ion torus proposal of Ref. [49]. 
In Section 4, we show the results of our simulations in Kerr and non-Kerr spacetimes and 
we study the correlation between the estimate of the spin and of possible deviations from 
the Kerr geometry. Summary and conclusions are reported in Section 5. Throughout the 
paper, we employ units in which Gn = c = I and the convention of a metric with signature 
(- + ++)• 
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2 Polish doughnut model in non-Kerr spacetimes 


Since we want to test the Kerr metric, we consider a backgronnd more general than the 
Kerr solution and that includes the Kerr solution as special case. In addition to the mass 
M and the spin parameter a* = a/M = where J is the BH spin angular momentum, 

the metric has one or more “deformation parameters”, which are used to quantify possible 
deviations from the Kerr geometry. The idea is to compute the spectrum of the ion torus 
in this more general background and see whether observations can potentially measure both 
the spin and the deformation parameter. The latter vanishes in the Kerr spacetime, but here 
it is assumed a free parameter to be inferred by observations. If the data require a vanishing 
deformation parameter, then the Kerr paradigm is conhrmed. As a non-Kerr background, 
we consider the Johannsen-Psaltis metric [53]. The simplest non-Kerr model has only one 
non-vanishing deformation parameter 63 . In Boyer-Lindquist coordinates, the line element 
reads 


ds^ = -\1- 


2Mr\ , X , 2 4aMr sin^ 0 _ , x , , , 

—\ h) dt -- [1 h) dtd(f) -|- 


S (1 + h) 

A -|- a^/isin^ 6 


dr^ + 




9 n 2a^Mrsm‘^ 9 
r^ + a^ + --r- 


2 „ a^(S-k 2Mr) sin^ 0 , 
sm"' 9 H-^^- h 


dcj)"^ 


( 2 . 1 ) 


where S = cos^ 9, A = — 2Mr + o^, and 


e^M^r 


( 2 . 2 ) 


The compact object is more prolate (oblate) than a Kerr BH with the same spin for 63 > 0 
(e 3 < 0); when €3 = 0, we recover the Kerr solution. 

The Polish doughnut model describes non-self-gravitating disks when the gas pressure 
cannot be neglected [51, 52]. This makes the disk geometrically thick and the particles of the 
gas do not follow the geodesics of the spacetime. The model only requires that the space- 
time is stationary and axisymmetric, and therefore its extension to non-Kerr backgrounds is 
straightforward, but new phenomena may show up [54]. The accretion flow is modeled as a 
perfect fluid with purely azimuthal velocity, so its energy-momentum tensor and 4-velocity 
are 

= {p + P)uf^u‘' + = (y, 0, 0, , (2.3) 

where p and P are, respectively, the energy density and the pressure of the gas. In what 
follows, we assume a polytropic equation of state 


P = , 


(2.4) 


where K is the so-called polytropic constant and n is the polytropic index. The specific 
energy of the fluid element, —ut, its angular velocity, Q = /u^, and its angular momentum 

per unit energy, I = —u^/ut, are 


Ut = 


9t(p 9tt9(f>(f> 

+ ‘^l9tij) + ^‘^9tt 


l9tt + 9tif> 
^9t(f) T 500 


5t0 T ^94>4' 

gtt + ^9t(j> 


(2.5) 


For a stationary and axisymmetric flow in a stationary and axisymmetric spacetime, I is 
conserved [51]. In the rest of this paper, we consider the simplest model with I = const. [52]. 
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It is useful to introduce the potential W 


lV(r,0) = - In 


9tt “1“ 9cl)(f> 

{gtt + 


( 2 . 6 ) 


The “center” of the accretion disk corresponds to the maximum of P and of W, say W = Wc- 
The surface of the torus is at some equipotential surface W = Wg where P = 0. Accretion is 
possible only when the fluid surface has some cusp: like the cusp at the Li Lagrange point 
in a close binary, the accreting gas can fill out the Roche lobe and then be transferred to 
the compact object. A cusp only exists when the fluid angular momentum I is between a 
maximum and a minimum value, say Imax and /min- If I > /max, the angular momentum is too 
high and the fluid cannot approach the central object. If / < /min, the angular momentum 
is too low and there is no stationary accretion flow: all the gas is quickly swallowed by the 
BH. In the Kerr metric and in many other backgrounds, /max corresponds to the angular 
momentum per unit energy of the marginally bound equatorial circular orbit, while /min is 
the angular momentum per unit energy of the marginally stable equatorial circular orbit. In 
some Johannsen-Psaltis spacetimes, there is no marginally bound equatorial circular orbit, 
but it is still possible to define /max, see Ref. [54] for more details. We just note that, in those 
backgrounds in which there is no marginally bound equatorial circular orbit, the accretion 
disk has two cusps, one above and one below the equatorial plane of the system. Instead of 
/ and W, we can use the dimensionless parameters 


A 


/ /r 


/max /ri 


w(r, 0) 


IK(r, 9) - IT, 
Wc-Wg 


(2.7) 


A can range from 0 to 1 by definition. The same is true for cj in the region occupied by the 
gas. 


3 Plasma model 

The accretion structure around SgrA* seems to be radiatively inefficient and advection dom¬ 
inated [48]. Following Refs. [49, 50], we employ the Polish doughnut model (reviewed in 
the previous section) to describe the macroscopic structure of the accretion disk and a gas 
pressure dominated optically thin 2-temperature plasma model for the microscopic physics. 
In this section, we just report the results, while more details can be found in Refs. [49, 50]. 
The model has 8 parameters: 2 related to the spacetime metric, 1 parameter associated to 
the position of the observer, and 5 parameters for the accretion flow and its microphysics: 


background metric spin parameter a* 

deformation parameter es 

observer viewing angle i 

ion torus dimensionless specific angular momentum A 

magnetic to total pressure ratio (3 

polytropic index n 

energy density at the center 
electron temperature at the center Tc 


We note, however, that A depends also on the spacetime, because the numerical values of 
/max and /min depend on the background metric. 
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We consider a polytropic index n = 3/2, which is suitable for a non-relativistic gas with 
negligible radiation pressure. The pressure P is the sum of the magnetic pressure Pm and 
the gas pressure Pgas, which are supposed to have a constant ratio, so 

B'^ 

P = Pm + Pgas , Pm = ■ ( 3 - 1 ) 


B is the intensity of the magnetic field. The gas is described by a 2-temperature plasma, so 
the gas pressure is the sum of the ion and electron contributions. The pressure P, the energy 
density p, and the electron temperature T can be written as [49, 50] 


where 


P = 
P = 

T = 


Kp^!^ , 




i^ = ^exp 
pc 


( W,-Ws 
V n + l 



(3.2) 

(3.3) 

(3.4) 


(3.5) 


With this set-up, we can compute the spectrum of an ion torus. We use the ray-tracing 
code described in Ref. [12]. The observed specific flux is given by 


Tobs(^obs) — / .^obs (^obs)dflobs — 19 .^e(^e)^^obs ; 


(3.6) 


where Jobs and Ig are, respectively, the specific intensity of the radiation as measured by the 
distant observer and in the rest frame of the accreting gas. dflobs = dXdY/d? is the element 
of the solid angle subtended by the image of the disk on the observer’s sky, X and Y are 
the Cartesian coordinates on the observer’s sky, and d is the distance of the BH from the 
observer. Jobs = follows from Liouville’s theorem, g is the redshift factor 


where is the 4-momentum of the photon, = (1,0,0,0) is the 4-velocity of the distant 
observer, and is the 4-velocity of the fluid element, k^/kt is a constant of motion along 
the photon path. The observer’s sky is divided into a number of small elements and the 
ray-tracing procedure provides the observed flux density from each element; summing up all 
the elements, we get T/bs- 

Contrary to the case discussed in Ref. [12], here the disk is geometrically thick and 
optically thin. The specihc intensity P is calculated by integrating along every photon path 
the emission and absorption contributions 




je(t'e) - a{Ve)P{ye 


d£, 


(3.8) 


where je and a are, respectively, the emission and absorption coefficient, d£ = u^kg,di is the 
infinitesimal proper length as measured in the rest frame of the emitter, and i is the affine 
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Figure 1. Left panel: spectrum of the ion torus of the default model {a/M = 0.5, €3 = 0, i = 60°, 
A = 0.3, /? = 0.1, n = 3/2, pc = 10“^^ g/cm^, Tc/T„ = 0.02). Right panel: contribution from 
every radiative mechanism to the default model, namely bremsstrahlung emission (blue solid line), 
Comptonisation of bremsstrahlung emission (violet dashed line), synchrotron radiation (red dotted 
line), and Comptonisation of synchrotron radiation (green dashed-dotted line). See the text for more 
details. 


parameter of the photon trajectory. In what follows we assume a = 0. The electromagnetic 
spectrum of the ion torus is produced by bremsstrahlung, synchrotron processes, and inverse 
Compton scattering of both bremsstrahlung and synchrotron photons off free electrons in the 
medium. Creation and annihilation of electron-positron pairs are neglected. The details can 
be found in [49, 50] and references therein. The emission coefficient is thus the sum of four 
contributions 


■ (3.9) 

The expression for these contributions are reported in [49, 50]. 

4 Simulations 

Fig. 1 show the spectrum of our default model. The background is described by the Kerr 
metric (es = 0) with spin parameter a* = 0.5. The viewing angle is i = 60°. Concern¬ 
ing the ion torus, we have assumed A = 0.3, /3 = 0.1, n = 3/2, p 2 = 10“^^ g/cvo?, and 
Tc/Ty = 0.02 where Ty is the virial temperature. The left panel shows the total spectrum, 
while the right panel shows the contributions from the emission mechanisms under con¬ 
sideration (bremsstrahlung, synchrotron emission, and inverse Compton scattering of both 
bremsstrahlung and synchrotron photons). 

The impact of the parameters of the model on the spectrum of the ion torus have 
been already discussed and illustrated in Ref. [49]. Here we just show the effect of the spin 
parameter a* and of our new parameter 63 describing possible deviations from the Kerr 
geometry. Figs. 2 and 3 show some spectra for different values of a* and 63 , respectively in 
the case A = 0.3 and 0.6. As in the case of the thermal spectrum of thin disks (see e.g. Fig. 2 
in Ref. [ 6 ]), the impact of a* and £3 looks similar and the spectrum could be degenerate with 
respect to these two parameters. However, the spectrum of a ion torus is more complicated 
and the conclusion is not obvious. 
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Figure 2. Effect of the spin parameter (left panel) and of the deformation parameter (right panel) 
on the spectrum of an ion torus. The other parameters are set to the value of the default model. See 
the text for more details. 


M 

I 10-'"' 


? -^-16 
^ MJ 


10 ' 


10 '“ 5 
10 ® 



_a/M=0.2 
a/M=0.5 
a/M=0.7 
a/M=0.9 


10 '-" 10 ' 
iz[Hz] 


10^ 



Figure 3. As in Fig. 2 for A = 0.6. 


To be a little bit more quantitative, we consider a reference model and we compare its 
spectrum with those produced in other models by varying the spin parameter, the deformation 
parameter, and possibly the specific fluid angular momentum A. We use the following function 
as an estimator 

5(a*,e3, A) = ^ 

i 


logT)(a^,e 3 , A) - logF- 
C log 


ref 


n 2 


(4.1) 


where is the observed intensity flux of the reference model at the frequency i, T)(a*, 63 , A) 
is the observed intensity flux of the model with (a*, 63 , A), and C is a constant that we have 
arbitrarily set to 0.1. All the other parameters are fixed to the value of the default model 
reported above. 

In Fig. 4, we show some contour levels of S in the case the parameters of the reference 
model are set to the values of the default model, whose spectrum is shown in Fig. 1. In the 
left panel, A = 0.3 is fixed. In the right panel, A is a free parameter and S is minimized with 
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Figure 4. Contour map of the function S in which the reference model has the parameters of the 
default model: = 0.5, = 0 (Kerr BH), = 0.3, = 0.1, = 3/2, g/cm^, 

and r/®^/T„ = 0.02. In the left panel, A = 0.3 is fixed, while in the right panel A is free and S is 
minimized with respect to it. See the text for more details. 



Figure 5. As in Fig. 4 with the reference model of a* = 0.9 and £3 = 0. In the left panel, A = 0.3 is 
fixed, while in the right panel A is free and S is minimized with respect to it. See the text for more 
details. 

respect to it, as one should do if A is unknown and to be determined by the fit. These plots 
clearly show the strong degeneracy between a* and 63 : it is impossible to get an estimate 
of one of these parameters without knowing the other from an independent measurement. 
Even if under the assumption of the Kerr background we could obtain a good constraint on 
the spin like 0.45 < a* < 0.55 (the contour line S = 2), actually we cannot exclude large 
deviations from the Kerr background and a possibly very different value of the spin a*. The 
impact of A is modest: when even this parameter is free, the constraints on a* and 63 are 
roughly the same as in the case of a hxed A. The degeneracy of the spectrum with respect 
to a* and 63 is very similar to that found in the measurements of the thermal spectrum of 
thin disks, see e.g. Ref. [15]. 

Fig. 5 is like Fig. 4, but now the reference BH has a* = 0.9 (all the other parameters 
are unchanged). A is fixed in the left panel, while it is free in the right panel. As with other 
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approaches, fast-rotating BHs are more suitable to test the Kerr metric and the constraints 
are stronger. However, the parameters a* and es are still very correlated and it is impossible 
to rule out large deviations from the Kerr geometry. 

5 Concluding remarks 

SgrA* is the supermassive BH candidate at the center of our Galaxy and a special source to 
test the Kerr metric. Unlike most of the other BH candidates, SgrA* may be soon studied 
with many different approaches and the combination of these measurements could provide 
an independent estimate of both the spin parameter and possible deviations from the Kerr 
solution. The general problem to test the Kerr metric is indeed that observational features 
of Kerr BHs may be reproduced by non-Kerr BHs with a different spin. If different measure¬ 
ments are sensitive to different relativistic effects, there is the possibility of breaking such a 
degeneracy and unambiguously testing the Kerr paradigm. 

In this paper, we have explored the possibility of testing SgrA* with the spectrum of 
its accretion structure. Following Ref. [49], we have employed the Polish doughnut model 
to describe an optically thin ion torus radiating via bremsstrahlung, synchrotron processes, 
and inverse Compton scattering. The model has 8 free parameters: 2 parameters related to 
the spacetime metric (spin parameter a*, deformation parameter 63 ), 1 parameter associated 
to the observer (viewing angle i), and 5 parameters to describe the ion torus (dimensionless 
specific angular momentum A, magnetic to total pressure ratio (3, polytropic index n, energy 
density at the center pc, and electron temperature at the center Tc). 

With the use of the ray-tracing code of Ref. [12], we have computed the spectrum of 
an ion torus with specific values of these 8 parameters. In the presence of good data of 
the spectrum of the accretion structure of SgrA*, one could potentially fit the observations 
and constrain the parameters of the model. In this explorative work, we have limited our 
attention to the spin, the deformation parameter, and to the fluid angular momentum, but 
a more accurate study should also check the correlation with the other parameters of the 
model. As in the case of the thermal spectrum of thin accretion disks, the spectrum of an 
ion torus is degenerate with respect to a* and 63 and it is impossible to determine one of the 
two parameters without an independent measurement of the other. In this sense, the sole 
observation of the spectrum of an ion torus cannot test the Kerr metric, even in the presence 
of high quality data and with the systematics under control. However, in the case of SgrA* 
it will be hopefully possible to combine several measurements. The approach studied in this 
paper may still provide interesting results if its constraints could be combined with those 
from other observations, and in particular from the observation of the BH shadow and from 
the detection of a pulsar orbiting with a period of a few months. 
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